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Abstract 



In the heavy quark hmit inclusive production rate of a heavy meson can be 
factorized, in which the nonperturbative effect related to the heavy meson can 
be characterized by matrix elements defined in the heavy quark effective theory. 
Using this factorization, predictions for the full spin density matrix of a spin-1 and 
spin-2 meson can be obtained and they are characterized only by one coefficient 
representing the nonperturbative effect. Predictions for spin-1 heavy meson are 
compared with experiment performed at e~^e~ colliders in the energy range from 
■y/s = lO.SGeV to ^/s = 91GeV, a complete agreement is found for D*- and B*- 
meson. For D** meson, our prediction suffers a large correction, as indicated 
by experimental data. There exists another approach by taking heavy mesons 
as bound systems, in which the total angular momentum of the light degrees of 
freedom is ^ and | for spin-1 and spin-2 meson respectively, then the diagonal parts 
of spin density matrices can be obtained. However, there are distinct differences 
in the predictions from the two approaches and they are discussed in detail. 

PACS numbers: 13.85.Ni,13.88.+e, 14.40.Lb,12.40.Nd 
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1. Introduction 

In the heavy quark hmit the momentum of a heavy meson containing one heavy quark can be 
approximated as the momentum of the heavy quark. This fact allows us to study properties of a 
heavy meson by starting from QCD directly. Activities of the study have leaded to the birth of the 
heavy quark effective theory(IIQET) |]T|,0], with this effective theory a consistent and systematic 
expansion in the inverse of the heavy quark mass can be preformed. HQET has been widely used 
in studies of decays of heavy hadrons, in this work we apply HQET for inclusive productions of a 
polarized heavy meson, where the heavy meson is a spin-1 or spin-2 particle. 

In its inclusive production a heavy meson is formed with a heavy quark Q and with other 
light degrees of freedom, the light degrees can be a system of light quarks and gluons. Because 
its large mass uiq the heavy quark is produced by interactions at short distance. Therefore 
the production can be studied with perturbative QCD. The heavy quark, once produced, will 
combine light degrees of freedom to form a hadron, the formation is a long-distance process, in 
which momentum transfers are small, hence the formed hadron will carry the most momentum 
of the heavy quark. The above discussion implies the production rate can be factorized, in which 
the perturbative part is for the production of a heavy quark, while the nonperturbative part is for 
the formation. For the nonperturbative part a expansion in the inverse of mq can systematically 
be performed in the framework of HQET. This type of the factorization was firstly used in parton 
fragmentation into a heavy hadron |0. 

In this work we will use the factorization to predict the spin density matrix of a spin-1- or spin-2 
heavy meson in its inclusive production. It turns out that the matrix is determined by perturbative 
coefficients and by one nonperturbative parameter, which is a ratio of matrix elements defined in 
HQET. The predictions are compared with experiment performed at e^e~ colliers in the energy 
range from ^/s = lO.SGeV to y/s = 91GeV, a good agreement is found for D*- and i?*-meson. 

The diagonal part of spin density matrices have been studied before [Q, in which the total 
angular momentum j of the light degrees of freedom in the heavy meson is taken as | and | for 
the spin-1- and spin-2 meson, respectively, i.e., j takes the smallest of possible values. It should 
be noted that this property of a heavy meson has been not derived from QCD, because solutions 
of eigenvalue problem of QCD hamiltonian in the heavy quark limit are not known. This property 
may be argued by that the system of light degrees of freedom with larger value of j is produced 
with a suppressed probability. For B* this property is experimentally supported by observing the 
ratio of the production rates of spin-O-and spin-1 states, while for charm mesons the ratio has 
a 20% deviation from the prediction based on the property. In our approach we will not take 
the property into account. In other word our results are derived on a more general ground and 
corrections from higher orders in tjiq^ can be systematically added. It is interesting to note that 
the obtained predictions of the spin density matrix of D* also agree with experiment well, in which 
one expects large corrections form higher orders of m'^. For a spin-2 heavy meson a parameter is 
introduced for its production in the approach in Q , and this is the only parameter characterizing 
the angular distribution of its decay. In our approach the corresponding distribution is predicted 
as an isotropic one. Experimentally this distribution is measured with a limited data sample for 
D**. Although the distribution is not well determined in experiment, the measured distribution 
does not agree with our prediction, indicating that the correction to our prediction made in the 
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heavy quark limit is substantial, this is also theoretically expected, because the binding energy of 
D** is quite large. We will discuss in detail the difference between the two approaches. 

In this work we consider the inclusive production at e^e~ colliders, where the initial beams are 
unpolarized. Our results can be generalized to the case with polarized beams, and also to inclusive 
production at hadron colliders or at ep-coUiders. It is straightforward to obtain predictions for 
other inclusive predictions. Our work is organized as the following: In Sect. 2 we give our result for 
spin-1 meson and compare our predictions with experiment. The differences between our approach 
and that in are discussed in detail. In Sect. 3 we present results for spin-2 meson. Sect. 4 is our 
summary. Throughout this work we take nonrelativistic normalization for heavy quarks and for 
heavy mesons. 



2. The density matrix for spin-1 heavy meson 



We denote a spin-1 heavy meson as H* , which contains one heavy quark Q. We study the 
inclusive production 

e+(p)+e-(-p)^iJ*(k)+X, (1) 

where the three momenta are given in the brackets. In the process we assume that the initial 
beams are unpolarized. We denote the helicity of H* as A and A = —1,0,1. All information 
about the polarization of H* is contained in a spin density matrix, which may be unnormalized or 
normalized, we will call them unnormalized or normalized spin density matrix, respectively. The 
unnormalized spin density matrix can be defined as 

R{\ A', p, k) = Y.{H*{\)X\T\e+e-) ■ (iJ*(A')X|T|e+e-)*, (2) 

X 

where the conservation of the total energy-momentum and the spin average of the initial state is 
implied. T is the transition operator. The cross-section with a given helicity A is given by: 

From Eq.(2) the normalized spin density matrix is defined by 

i?(A,A',p,k) 

^"''P-''> = E.fi(A.A.P.k) - W 

It should be noted that the normalized spin density matrix is measured in experiment. 

The contribution to -R(A, A', p, k) can be characterized by Feynman diagrams. With the fac- 
torization discussed in the introduction, each diagram can be divided into two parts. For the 
process in Eq.(l) the contribution at tree-level is given in Fig.l, where we divide the diagram 
with a horizontal broken line into two parts. The upper part with the black box represents the 
formation, while the lower part represents the inclusive production of a heavy quark Q, which is 
not necessarily on-shell here. At loop-level it is possible that there are gluon exchanges between 
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quark lines and between the lower part and the black box. The effect of the exchanges between 
the lower part and the black box is at higher orders in w^g^, we will neglect it in this work. The 
exchanges of gluons between quark lines will in general bring some infrared singularities in the 
lower part. We can show that at the leading order of m'^ the heavy quark lines connecting the 
lower- and upper part represent on-shell quarks. Hence, the lower part with possible gluon ex- 
changes is responsible for inclusive production of an on-shell quark Q, and it is free from infrared 
singularities. 



q 



7, 



7, 



Fig.l: The diagram for the contribution at tree-level. The line is for heavy quarks. The vertical 
broken line is the Cutkosky cut. 



With the diagram the density matrix can be written as 
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i?(A,A',p,k) = J ^T,,(g)r,,(A,A',g,fc), (5) 

where Tji{q) is for the lower part and can be calculated with perturbative theory, Tij{\, A', q, k) 
represents the formation of the heavy meson, which is defined as: 

r,,(A,A',g,^) = / d^xe-^^-''^(0|g,(0)|i/*(A')X)(M*(A)|g,(x)|0). (6) 
■' X 

The indices i.j stand for Dirac- and color indices, they connect the lower- and the upper part 
through the vertical quark lines. Q{x) is the Dirac field for the heavy quark Q. In the heavy quark 
limit the heavy meson carries the most momentum of the heavy quark, the difference between the 
momentum of the heavy quark and that of the heavy meson is at order of A.qcd- This scale Aqcd 
is much smaller than rriQ and ^/s, which are scales appearing in the lower part of Fig.l. This 
fact indicates that we can expand A', q, k) in Aqcd- With this expansion for Tij{\, A', q, k) 

it results in that the density matrix given in Eq.(6) is expanded in rriQ^. For doing this we use 
HQET, in which the Dirac field is expanded in the field of HQET: 

QW = e-"0«-{i + J^n.flT)/..M + O(^) 

+ (terms for anti — quark), (7) 

where v is the velocity of H*, h^{x) is the heavy quark field in HQET, Df^ = — v ■ Dv^ and 
is the covariant derivative. Substituting this into rjj(A, A', g. A;), the integrand depends on x 
through hv{x) and through a exponential factor, the x-dependence through hv{x) is controlled by 
the scale at order of Aqcd, at leading order this dependence can be neglected. We then have: 

ri^{X,X',q,k)^{2n)^S\q-mQv)Y,{0\{K{Q))i\H*{X')X){XH*{^^^^^ (8) 

X 

where the ■ ■ ■ stand for neglected terms, whose contributions to the density matrix are suppressed 
by powers of ttiq^. With a detailed examination one can show these contributions are of order of 
itiq^ or of higher orders. At the leading order the momentum of H* is approximated by mQV. To 
analyze the structure of Fij in the above equation it is convenient to work with the rest frame of 
H*. which is obtained from the moving frame only by a Lorentz boost. We choose the 2;-direction 
to be the moving direction of H* and denote the polarization vector of H* in the rest-frame as 
e(A). In this frame we can define a creation operator for H*: 

|//*(A)) = at(A)|0) = e(A)-at|0). (9) 

In the rest frame the field has two non zero components and can be written as: 



hv{x)=[^^yj. (10) 
With this notation the matrix in Eq.(8) takes the form: 
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Ex mK{0))i\H*{\')X){XH*{X)\ChM)M 

f l + l-v / (0|^at(A')a(A)V^t|o), o\ 1 + 7-^ ^ 
" 2 ■ 0, J ■ 2 



(11) 



The matrix element {0\ipa^{X')a{X)ilj^O) is diagonal in color space, the matrix structure in spin 
space can be decomposed into the unit matrix and the Pauli matrices cr. For the decomposition 
we define two operators: 

0{H*) = lTii;alaiij\ Os{H*) = ^1iaii,a]aki^hiju, (12) 
1/ 

where Sij^ is the totally antisymmetric tensor. With these operators the matrix can be written: 
(0|V'at(A')a(A)^t|o) = i(0|O(if*)|0)6*(A) ■ 6(A') -I + l^- [e\\) x e(A')] ■ (0|O.(if*)|0), (13) 

where / is a 2 x 2 unit matrix in the spin space. It should be noted that the factor (1 + 7 ■ f )/2 
in Eq.(ll) is a projector and can be written as Es w(j9, s)'u(p, s), where u{p,s) is the spinor of 
the heavy quark Q in the rest frame. This means that the projector projects the quark state 
represented by the vertical quark lines in Fig.l to the on-shell state. Therefore, the lower part is 
just for inclusive production of a on-shell heavy quark Q. The two matrix elements {0\O{H*)\0) 
and {0\Os{H*)\0) are universal in the sense that they do not depend on how Q is produced, they 
describe the long-distance process of the formation of H*. Using Eq.(8), Eq.(ll) and Eq.(13) we 
can write the result for the unnormalized spin density matrix as: 

P(A, A',p,k) = la(p,k)(0|O(if*)|0)6*(A) • 6(A0 + ^b(p,k) ■ [e*(A) x e(A')] • {0\O,{H*)\0). (14) 

The quantities a(p, k) and b(p, k) characterize the spin density matrix of the heavy quark Q 
produced in the inclusive process: 

e+(p) + e-(-p)^Q(k,s)+X (15) 

where s is the spin vector of Q in its rest frame and the rest frame is related to the moving frame 
only by a Lorentz boost. The unnormalized spin density matrix Rq{s, p, k) of Q can be defined by 
replacing H*{X) with Q(k, s) in Eq.(2). This matrix can be calculated with perturbative theory 
because of the heavy mass. The result in general takes the form 

^Q(s,P,k) =a(p,k) +b(p,k) -s (16) 

where a(p, k) and b(p, k) are the same in Eq.(14). The results at tree-level for a(p, k) and b(p, k) 
may be found in 0. With Eq.(14) we obtain the normalized spin density matrix: 



P(P, k) = ^ 



(I + P3, 

-P-, I, -P+ I, (17) 
^ 0, 
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with 



&3(p,k) {0\OsiH*)\0) 
a(p,k) ■ {0\OiH*)\0) 



fei(p,k)±z62(p,k) (0|O,(if*)|0) 
V2a(p,k) ■ (0|O(i/*)|0)' 



(18) 



The indices of matrices in Eq.(17) run from -1 to 1. 

Before confronting experimental results we give the following comments to our results: 

1) . If one uses perturbative theory to calculate the coefficients in Eq.(14) or Eq.(16), one 
will encounter terms with large logarithm \\i{mQ/^). One needs to resume these terms. One 
way to resume is to use the well known factorization formula for the process in Eq.(l), where 
the heavy meson H* is produced through parton fragmentation. The definitions of the spin- 
dependent fragmentation functions can be found in Then one applies the approach here 
to the fragmentation functions as used in |^. These functions can be calculated at the energy 
scale /i = mq with perturbative QCD. Through renormalization group equations one obtains the 
functions and then the coefficients at the energy scale /i = ^i. Before the resummation the 
coefficients at tree-level are singular in the energy of H* , but the ratio P3 and P± are regular. 
After the resummation is performed, the coefficients become regular. 

2) . If the parity is conserved and T-odd effects can be neglected, which are proportional to 
mg/v^, then b(p, k) = 0, i.e., the heavy quark is unpolarized, it leads to that H* is unpolarized 
too. We obtain in this case poo = P-1-1 = Pii = 1/3, and all nondiagonal elements of the spin 
density matrix are zero. It should be noted that in general the nondiagonal elements can be 
nonzero even if the parity is conserved, e.g., a light vector meson can have tensor-polarization, 
which partly corresponds to the matrix element p_aa ^01 A = ±1. Unfortunately, the spin of H* is 
analyzed in experiment with its parity-conserving two-body decay, where the polarization of the 
decay products is not observed. This leads to that the polarization of the heavy quark, i.e., the 
contributions from b(p, k) to the spin matrix element, can not be determined in experiment, this 
is the so-called "no-win" theorem 0]. 

3) . Our approach can be used not only for the process in Eq.(l), but also for other processes. 
For other inclusive productions the same analysis can be done, where one needs to replace the 
lower part in Fig.l and Tji{q) in Eq.(5) with those for the corresponding production of Q, while 
the upper part in Fig.l and Tij in Eq.(5) remain the same. With the analysis presented here it 
is straightforward to obtain predictions for other inclusive productions, and predictions always 
takes the form as given in Eq.(14) or in Eq.(17), where the coefficients a(p, k) and b(p, k) should 
be replaced by those characterizing the spin density matrix of Q in the corresponding processes. 
Without knowing these coefficients we can always conclude that in any inclusive production of H* 
we have poo = 1/3 and p_ii = pi_i = in the heavy quark limit. These predictions may be tested 
in other processes, e.g., in inclusive production at electron-proton colliders or at hadron-hadron 
colliders. It is to the author's knowledge that the spin-measurements for B*, D* and D** have 
been done only at e"'"e~-colliders by those experimental groups, whose results are listed in this 
work. The list may be incomplete. It would be interesting to test our predictions in experiment 
at an electron-proton- or a hadron-hadron collider. 

The experiments to measure the polarization of B* are performed at LEP with ^/s = Mz by 
different experiments groups. To measure the polarization the dominant decay B* 7!? is used, 
where the polarization of the photon is not observed. Because the parity is conserved and the 
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distribution of the angle between the moving directions of 7 and of B* is measured, one can only 
determine the matrix element poo- If we denote 6 is the angle between the moving directions of 
B* and of 7 in the i?*-rest frame and is the azimuthal angle of 7, then the angular distribution 
is given by Ws'^B-yid, (j)) oc J2\\' Pwi^w - Yix{9, (p)Y^x'id^ 'P))- Integrating over and using 
our result poo = 1/3, the distribution of 9 is isotropic. In experiment one indeed finds that the 
distribution is isotropic in 6. The experimental results at ^/s = Mz are: 

Poo = 0.32 ± 0.04 ± 0.03, DELPHI [8] 
Poo = 0.33 ±0.06 ±0.05, ALEPH 0, 

Poo = 0.36 ±0.06 ±0.07, OPAL (19) 

These results agree well with our prediction poo = 1/3. It is interesting to note that the spin 
alignment is also studied with the LUND string fragmentation model |jll| implemented with 



JETSET in |T^, and the result poo = 0.567 is obtained, which is not in agreement with 
experiment. In 0] the B*- and B meson is taken as a bound state of a 6-quark and other light 
degrees of freedom, and these light degrees have the total angular momentum j = \- Because the 
parity is conserved in the formation of B* , the light degrees are produced with equal probability 
for ja = ±|. With this argument one obtains the probability for a left-handed 6-quark to form a 
B*- and B meson is: 

P{B\\ = -1)) : P(5*(A = 0)) : P{B*{X = 1)) : P(5) = ^ : ^ : : ^. (20) 

One then also obtains poo = 1/3 in agreement with experiment and with our results. One can also 
obtains the ratio of the production rates 



V 



(j{B*) 3 



a(P*)±a(P) A'' 
this prediction is well tested with experiment: 



(21) 



Py = 0.72 ±0.03 ±0.06, DELPHI [8] 

Pv = 0.771 ± 0.026 ± 0.07, ALEPH [| . (22) 

However the prediction of the ratio for c-fiavored mesons has a deviation at 20% level. It should 
be noted that in only the diagonal part of the density matrix is predicted, while in our work 
the predictions are given for the complete matrix. 

The polarization measurement for Z}*-meson has been done with different ^/s, in some exper- 
iments the non-diagonal part of the spin density matrix has also been measured by measuring 
azimuthal angular distribution in D* decay, where the decay mode into two pseudo-scalars, i.e., 
D* — > Dn, is used. Denoting 6 as the angle between the moving directions of D* and of vr in the 
D*-Test frame and (p as the azimuthal angle of vr, then the angular distribution of tt is given by 
WD'^Dn{0,(l)) oc J2xx' P\x'yix{0,(p)Y{y{9,(l)). Integrating over and using our result poo = 1/3, 
the distribution of 6 is again isotropic. The experimental results are summarized in Table 1 and 



also partly summarized in |]I4 
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Table 1. Experimental Results for D* 



Collaboration 



^/s in GeV 



Results 



CLEO 15 



HRS |T6| 



TPC |T7| 



SLD ]TB| 
OPAL 



10 



10.5 
29 



29 



91 
91 



Poo = 
Poo = 

pl-1 

PlO = 

Poo = 

Pl-1 

PlO = 

Poo = 
Poo = 

Pl-1 



0.327 ±0.006 

0.371 ±0.016 
= 0.04 ±0.03 

0.00 ±0.01 

0.301 ±0.042 ±0.007 
= 0.01 ±0.03 ±0.00 

0.03 ±0.03 ±0.00 

0.34 ±0.08 ±0.13 

0.40 ±0.02 ±0.01 
= -0.039 ±0.014 



From Table 1. we can see that the poo nieasured by all experimental groups is close to the 
prediction poo = |; the most precise result is obtained by CLEO, its deviation from the prediction 
is 2%, the largest deviation of the prediction is from the result made by OPAL at y/s = 90GeV, 
it is 20%. In general, poo depends on the energy of H*. Our results in Eq.(17) give that poo is a 
constant in the heavy quark limit, or the energy dependence is suppressed by rriQ^. In experiment 
only a very weak energy dependence is observed, e.g., in CLEO results [^]. From our results pi_i 
is exactly zero in the heavy quark limit, the results from TPC and from HRS are in consistent 
with our result, a non zero value is obtained by OPAL, which has a 3a deviation from zero. These 
deviations may be explained with effects of higher orders in m~^, these effects are expected to 
be substantial, because rric is not so large. It is interesting to note only results from OPAL at 
■s/s = 91GeV have the largest deviations from our predictions, while results from other groups 
agree well with our predictions. At ^/s = 10.5GeV or 29GeV, the effect of the Z-boson exchange 
can be neglected, hence the parity is conserved. We obtain pio = 0. This prediction is also in 
agreement with the experimental result made by TPC and by HRS. 

Since our results are derived without knowing the total angular momentum j of the light 
degrees of freedom in the heavy meson, the agreement of our results with experiment can not be 
used to extract the information about j from the experimental data in Eq.(19) and in Table 1., 
although Poo = 1/3 can also be obtained by taking j = 1/2. One way to extract j may be to 
measure the difference between pn — p_i_i, but it seems not possible, because the polarization 
of H* is measured through its parity-conserved decay and the polarization of decay products is 
not observed in experiment. In the heavy quark limit, the nondiagonal element pi_i and p_ii are 
zero, while the other nondiagonal matrix elements are nonzero if the parity is not conserved and 
the initial state is unpolarized. At higher orders in mq^ this can be changed, e.g., H* can have 
tensor polarization. 



3. The spin density matrix for spin-2 heavy meson 



In this section we consider a spin-2 heavy meson H**, which contains one heavy quark Q. We 
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study the inclusive production 



(P) + e" 



H**{k)+X. 



(23) 



With the detailed analysis in the last section it is straightforward to obtain the spin density 
matrix. To avoid introducing too many notations, we will use some notations used in the last 
section. The used notations in this section are referred to H**. We denote the unnormalized spin 
density matrix i?(A, A', p, k), which is defined by replacing H* with H** in Eq.(2). The hehcity A 
takes the value from -2 to 2. We work in the H**-Test frame, which is related to the moving frame 
only through a Lorentz boost. The 2;-direction is chosen as the moving direction of H**. In the 
rest frame the polarization tensor of H** is ejj (A), it has the following properties: 



eij{X) ^ eji{X), eii(A) = 0, 

A ^ d 

We introduce a creation operator for H** in its rest-frame: 

\H*%X)) = e,,{X)al\0), 



(24) 



(25) 



where aj is symmetric and trace-less. In analogy to Eq.(12) for the spin-1 case two operators can 



be defined: 

0{H**) = ^TVV'a^a,,•V'^ 0,{H**) = ^TVa,V'sV^fcV'^£i,fe- (26) 
With these operators the matrix for H** corresponding to that for H* in Eq. (13) can be written: 



E 

X 



H**{X')X){XH**{X)\^^;^0) = -{0\O{H**)\0)e*^{X)e,j{X') ■ I 



+-{OmH**)\0)e*,{X)eM{X')eiji ■ ai, 



(27) 



the matrix is diagonal in color space. Predictions can be obtained by replacing the matrix element 
in r.h.s. of Eq.(ll) with that given above and perform the same calculation as that in the last 
section. The unnormalized spin density matrix for H** then reads: 

R{X,X',p,k) = ia(p,k)(0|O(i/**)|0)6^(A)6,,(A') + ^6,(p, k)(0|O,(i/**)|0)e*,(A)6«(A')%;, (28) 

where the coefficients a(p, k) and 6j(p, k) for i = 1, 2, 3 are the same as those defined in Eq.(16). 
With this matrix we obtain the normalized spin density matrix for H**: 



pip, k) 



/ 1 + ^3, 


1 p 


0, 


0, 


\ 


0, 


1 + i^3, 




0, 





V3p 


1, 


2"-^+' 





0, 


0, 


a/3 p 
2 


1 - 1^3, 


- 1 p 


V 0' 


0, 


0, 


Lp 





(29) 
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with 



^ &3(p,k) (0|O.(g-)|0) ^ b^{p,k)±^b2{p,k) {0\O,{H**)\0) 

' a(p,k) (0|O(iJ-)|0)' ^ v^a(p,k) (0|O(/f-)|0) ' ^ ^ 

The indices of matrices in Eq.(29) run from -2 to 2. The results are similar to those for the spin-1 
case. The correction to the predictions in Eq.(29) is of order of mg^ or of higher orders. The same 
comments made in the last section for spin-1 heavy meson also apply here. 

It would be interesting to compare the predictions with experiment. Unfortunately, there is 
no experimental data for B** meson, there is a candidate -B}(5732) observed in experiment, but 
its spin is not determined. For c-flavored meson, a spin-2 meson noted as 1^2(2460) is observed, 
and the information of its spin is analyzed by ARGUS |[TP| with a limited data sample. Through 
the angular distribution of the decay D2 —>■ Dtt, it is found that there is no significant population 
of |A| = 2 states. For the decay there is one form-factor in the decay amplitude, the angular 
distribution can then be written as: 

(d; ^ Dn) ~ 5:pAA'r2A(^,0)r;v(^,0), (si) 



Td(j)d cos 6 



XX' 



where 6 is the angle between the moving directions of D2 and of vr in the D2-rest frame, (p is the 
azimuthal angle of vr. With our results we obtain: 

(D;^Dn) = l, (32) 



TdcosO' ^ ' 2 

i.e., the distribution is isotropic. In the approach in the distribution is predicted as: 

— -{D; Dtt) = -[1 + Scos^^ - 6cu3/2(3 cos^ ^ - 1)], (33) 

i a cos 4 

where an unknown parameter is introduced, this parameter is interpreted as the probability 
to produce a system of light degrees of freedom with the total angular momentum j = 3/2 and 
_73 = 3/2, and this system is combined with a heavy quark Q to form H**. ARGUS has found 
that the data favored a small a;3/2, indicating that a small portion of |A| = 2 states is produced. 
A re-analysis of the data is given in [|]. Assuming 0^3/2 > 0, it is found [Q: 

uj-i/2 < 0.24, 90% C.L. (34) 



A model estimation also gives a small UJ3/2 pOf. With our prediction we obtain that the distri- 



bution is isotropic. This leads to UJ3/2 = 1/2. It seems that experimental results are against our 
predictions. However, one should note that the correction from higher orders of can be more 
significant than that in the case of D*, because the excited state D2 is heavier than D*. The 
effects of higher orders may be estimated by the mass difference: 

MjP*) - MjDo) 
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this indicates that the large correction of our prediction (^3/2 = 1/2 exists. An extension of our 
present work to the next-to- leading order in is needed to analyze this effect, including effects 
from higher orders the distribution will be not isotropic. At leading order of the expansion in 
mg^ we obtain the relation pxx + P-\-\ = 2/5, this leads to that the distribution in Eq.(32) or 
(33) is isotropic. If we add corrections at higher orders of ttiq^, the relation will be changed, 
the coefficient 2/5 receives corrections which depend on A. This will change the distribution as 
predicted in Eq.(32) and it becomes non-isotropic. It needs to be study further. For B** meson 
the correction is expected to be small because rrih is large, hence the predictions are more accurate. 
It would be interesting to test our prediction in future experiment with B**. 

4. Summary 

In this work we have studied the polarization of a heavy meson in its inclusive production, 
where the heavy meson contains one heavy quark. A QCD factorization is performed for the 
spin density matrix of a spin-l-and spin-2 meson. In the heavy quark limit, the spin density 
matrix takes a factorized form, in which the nonperturbative effect related to the heavy meson is 
represented by two matrix elements defined in HQET, while the perturbative part is related to 
the spin density matrix of the heavy quark in its inclusive production. Hence, a detailed relation 
between the spin density matrices of the heavy meson and of the heavy quark is established. 

Without the knowledge of the two matrix elements defined in HQET, some elements of the 
spin density matrix can already be predicted. These elements are well measured in experiment. 
We find that the measured element for i?*-meson at LEP with = 91GeV agrees well with 
our prediction. For D*-meson the polarization is measured at e^e~ colliders with -^i = 10.5, 29 
and QlGeV, the experimental results also agree with our predictions, except a small deviation 
with OPAL data, the deviation may be due to corrections from higher orders in m~^. In general 
one expects these corrections are large because rric is not so large, the agreement for D*-meson 
indicates that HQET can be well used to predict the polarization of D*-meson. For D** meson 
not many experimental data are available. The results of a measurement with large errors does 
not agree with our result. The disagreement may be explained with large corrections from higher 
orders in these large corrections for D** meson are more likely than those for D* meson, 

because the binding energy of D** is larger than that of D* meson. 

The polarization of a heavy meson in its inclusive production has also been studied in [^], 
where one takes the total angular momentum of light degrees of freedom in the heavy meson to 
be the lowest of the possible values. With this the diagonal part of the spin density matrix can be 
predicted. In our work we only use the QCD factorization concept by employing the expansion in 
the inverse of the heavy quark mass, and we obtain results not only for the diagonal part, but also 
for the nondiagonal part of the spin density matrix. For the spin-1 case we obtain the same matrix 
element poo; while for the spin-2 case, our results are different than those in it leads to that 
the prediction for the angular distribution of the decay D** Dn is different than that predicted 
in 0. However, large corrections to our prediction are likely possible, it will be interesting to 
test predictions for B** meson in future experiments, because corrections in the case with B** are 
expected to be small. 
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Although we have given in this work detailed predictions for inclusive production of a polarized 

heavy meson at a e+e" collider, our approach can be easily generalized to other inclusive produc- 
tions, testable predictions can be made without a detailed calculation, for example, in inclusive 
production of B* at an electron-hadron- or a hadron-hadron collider we always have the prediction 
Poo = 1/3 and = = in the heavy quark limit. 
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